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Method of extending hyperfine coherence times in Pr 3+ : Y 2 Si05 
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In this letter we present a method for increasing the coherence time of praseodymium hyperfine 
ground state transitions in Pr 3+ :Y2SiOs by the application of a specific external magnetic field. 
The magnitude and angle of the external field is applied such that the Zeeman splitting of a hy- 
perfine transition is at a critical point in three dimensions, making the first order Zeeman shift 
vanishingly small for the transition. This reduces the influence of the magnetic interactions between 
the praseodymium ions and the spins in the host lattice on the transition frequency. Using this 
method a phase memory time of 82ms was observed, a value two orders of magnitude greater than 
previously reported. It is shown that the residual dephasing is amenable quantum error correction. 

PACS numbers: 



There is growing interest in the use of nuclear spin 
states associated with dopant ions in a solid state host 
to store and manipulate quantum information^, 0, 0- 
These applications require the relevant spin transitions 
to have long coherence times. This can be challenging 
to achieve in a solid state system due to magnetic in- 
teractions between the dopant ion and spins within the 
host. A spin free host can be utilized to minimize these 
interactions 0, however not all dopant species of interest 
are chemically compatible with such hosts. An example 
of a class of dopants where no satisfactory spin free host 
has been identified, are the rare earth ions. Because of 
the potential to manipulate their spin states optically, 
the use of rare earth ions has been proposed in a num- 
ber of quantum information applications^, In this 
letter we investigate the decoherence in the ground state 
hyperfine transitions in Pr^ + \Y2SiO5 and demonstrate a 
new method for increasing their coherence times. This 
technique is expected to be applicable to a wide range of 
spin systems. 

p r 3 +:Y 2 Si05 was chosen to investigate since along 
with its use in a quantum computing architecture having 
been proposed by Ichimura et al.[l|, it has already been 
used in slow and stopped light demonstrations^,!^]. The 
longest coherence time (T2) stated in the literature for 
a ground state hyperfine transition in Pr 3+ '.Y^StO^ is 
500/is6j. This is likely to be insufficient for practical 
quantum computing applications given that the hyper- 
fine transition frequency is ^lOMHz, limiting the Rabi 
frequencies of any driving field to less than a few MHz so 
as to be transition specific. This only allows of the order 
of 1000 operations to be completed within T 2 . 

Yttrium orthosilicate has symmetry given by the C 2h 
space group with two formula units Y2SiOs per trans- 
lational unit. This gives four different sites at which the 
Pr can substitute for Y. The four sites can be divided 
into two pairs. The two pairs labelled 'site 1' and 'site 2' 
@ have different crystal field splittings and hence differ- 



ent optical and hyperfine frequencies. The members of 
each pair, labelled 'a' and 'b', are related to each other 
by the crystal's C 2 axis and their optical and hyperfine 
frequencies are degenerate in zero field. The Pr sites in 
Y2SiC>5 have electronic singlet ground states with hyper- 
fine ground state interactions described by the following 
Hamiltonian: 



H = BMI + IQI 



(1) 



where B is the magnetic field vector, I is the vec- 
tor of nuclear spin operators, M is the effective Zeeman 
tensor combining nuclear and electronic Zeeman interac- 
tions and Q is the effective quadrupole tensor combining 
the quadrupole and second order magnetic hyperfine in- 
teraction, known as the pseudo-quadrupolc 9|. The one 
naturally occurring isotope of Pr has a nuclear spin of 
5/2. M and Q for site 1 were recently determined for 
Pr 3+ :Y 2 Si0 5 to be: 



Q = 



M = 
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where E = 0.5624M7iz and D = 4A45QMhz, 
{g x ,g y ,gz) — (2.86, 3.05, 11.56)kHz/G, and the Euler an- 
gles are (a, £,7) = (-99.7,55.7,-40) 10]. These values 
are for the crystal aligned with the C2 axis in the y direc- 
tion, and the z axis is the direction of linear polarization 
of the praseodymium optical transitions. These tensors 
are highly anisotropic due to the low symmetry of the 
site. 

The dominant dephasing mechanism for the Pr hyper- 
fine ground states in Y2SiOs is due to magnetic inter- 
actions with the Y nuclear spins in the host. Yttrium 
has a nuclear spin of 1/2 and a gyromagnetic ratio of 
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FIG. 1: Two and three pulse echo sequences. 7r pulse duration 
was 40/is and the laser was turned off 150>s before the pulse 
sequence and turned on 100>s before the echo 
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FIG. 2: Spectrum of hyperfine transitions showing experi- 
mentally measured points as the field is increased toward the 
critical point 



g n = 2QQHz/G. Nearest neighbor Y ions induce a mag- 
netic field at the Pr site of the order of 0.1G. The direct 
magnetic dipole-dipole interaction between the Pr and 
a nearest neighbor Y ion is of the order of hundreds of 
Hertz. This interaction can induce decoherence via two 
mechanisms. The first mechanism is that the Y ions can 
exchange spin with each other, resulting in a fluctuat- 
ing magnetic field at the Pr site. For resonant Y nuclei 
the rate of exchange between two nearest neighbor ions 
is of the order of tens of Hertz. The second mechanism 
involves the excitation of nearly degenerate transitions 
where one or more Y spin flips are induced along with 
the desired change in the Pr hyperfine state. 

To suppress these dephasing mechanisms for a given 
Pr ground state hyperfine transition we propose apply- 
ing a specific magnetic field such that transition's linear 
Zeeman shift about this field value is zero. The main re- 
quirements for this technique are the existence of a zero 
field splitting in the spin states and the ability to apply 
a sufficiently strong magnetic field such that the Zeeman 
splitting is comparable to the zero field splittings. This 
method should be applicable to many I > 1 spin systems. 

Using equation [2] and the parameters for the site 1 
given above, all the magnetic field values resulting in a 
zero first Zeeman shift were identified. The case chosen 
to experimentally investigate was at a magnetic field of 
B C p = {732,173, -219} G on the m 7 = +1/2 <-> +3/2 
transition at 8.63 MHz as it was found to have the small- 
est second order Zeeman shift. Around this magnetic 
field value the transition energy as a function of mag- 
netic field has a turning point in the y and z axis while 
the x axis has a slow inflection point. As the required 
magnetic field was not along the crystal's C2 axis the de- 
generacy between site la and site lb was lifted with only 
the site la ions being at the critical point. 

Raman heterodyne two and three pulse spin echoes 
were used to investigate the decoherence and spectral dif- 
fusion in the ground state hyperfine transitions 6J (figure 
The experiment was performed using a Coherent 699 
frequency stabilized (1MHz FWHM) tunable dye laser 
tuned to the 3 H4 — 1 Z?2 transition at 605.977nm. The 
laser was gated using a 100MHz AOM such that there 
was no laser radiation applied to the sample during the 



RF pulse sequence (figure ^) . The laser power incident 
on the crystal was 40mW, focused to ~100/im and a RF 
Rabi frequency of Q RF = 300kHz. The Pr 3+ :Y 2 Si0 5 
(0.05% concentration) crystal was held at ~ 1.2K for 
the duration of the experiment. The laser prepared a 
population difference in the excited ions for 5s before the 
pulse sequence and was scanned over 1.2GHz of the opti- 
cal transition to avoid hole burning effects. The laser is 
off during the RF pulse sequence to minimize coherence 
loss from optical pumping. 

The magnetic fields were supplied by two supercon- 
ducting magnets one along the z axis the other in the hor- 
izontal plane. The sample was then rotated by 13 ± 0.5° 
to provide the correct ratio of fields along the x and y axes 
for the critical point in magnetic field space. Corse ad- 
justment of the field was done by comparing the Raman 
heterodyne spectrum to theoretical simulations as shown 
in figure|21 Fine adjustment of the magnetic field utilized 
perturbing coils along each axis and associated lock-in 
amplifiers to perform field sensitivity measurements on 
the frequency of the mi = +1/2 «-> +3/2 transition. 
The rotation of the sample and the current driving the 
magnets were iteratively adjusted to minimize the sen- 
sitivity of the desired transition. Final adjustments of 
the field values were made by minimizing the echo decay 
rate. The magnetic field values were accurate to within 
±0.1 G of the critical point. 

The two pulse echo data (figure |3J) shows three echo 
sequences at the critical point magnetic field configu- 
ration with the zero field echo sequence for reference. 
Since the field configuration is transition specific the 
mi = +1/2 <-> +3/2 transition at site la is at a criti- 
cal point while the two other transitions are not. 

With the application of the field for all the transitions 
there was observed a significant increase in the coher- 
ence time. For the mi = +3/2 <-> —3/2 transition T2 
was 5.86ms while for the mi = +1/2 <-» +3/2 transi- 
tion at site lb it was 9.98ms (figure EJ. Applying a field 
lifts the degeneracy of the Y spin states and inhibits Pr 
transition involving single Y spin flips[ll|. This can be 
understood by considering the Y quantization axis. At 
zero field the Y quantization axis is locally defined by the 
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FIG. 3: Two pulse echo sequences taken at the critical point 
magnetic field Bcp = {732, 172, —219} G for transition mi — 
+1/2 <-» +3/2 at site la (+). The mi = +1/2 <-» +3/2 site 
lb transition (x); the mj = —3/2 <-» +3/2 transition (o) and 
the zero field T2 = 500/^s decay line shown for comparison. 



Pr ion. Therefore a change in the spin state changes the 
quantization axis, mixing the Y spin states and result- 
ing in a high spin flip probability. With an applied field, 
large compared to the Pr-Y coupling, changes in the Pr 
spin state only weakly perturb the quantization axis. 

The echo decay for the mj = +1/2 <-> +3/2 at site la is 
significantly slower than for the other two transitions and 
is clearly not described by a simple exponential. Fitted 
to this decay is an exponential with a quadratic time de- 
pendence, equation[3] This coherence decay is indicative 
of a spectral diffusion process described by a Lorentz ker- 
nel with a width increasing linearly with time at a rate of 
47.3-ffz/s|l2j. The phase memory, Tm is used to charac- 
terise the decoherence rate of the system. In the present 
case the Tm was determined to be 82 ms. 



1(2*) = I exp 



2* 
Tm 



(3) 



To investigate the diffusion process over longer times 
scales a series of three pulse echo measurements were 
made. Within each of the series shown in FigQ]the delay 
between the first two pulses (tj.) is held constant while 
the delay between the second and third pulses (T2) is 
varied. 

A comparison is made with the same transitions as 
used in the two pulse echo study. Comparing echo se- 
quences with the same t\ on the different transitions 
shows a clear difference between the critical point and 
the other two transitions (figure^). The mj — —3/2 <-> 
+3/2 and the site lb mi = +1/2 <-> +3/2 transitions 
have a clearly biexponential decay: a fast initial decay, 
reducing to a slow decay for T2 > 200ms. With a t\ of 
2.5ms no spectral diffusion was observed on the critical 
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FIG. 4: Three pulse echo sequences taken at (a) mi = 
+ 1/2 <-> +3/2 transition at site la(+) and site lb(.) mi = 
—3/2 <-> +3/2 transition (o), both traces were taken with 
n — 2.5ms. The delay T2 is measured on the x axis, (b) Crit- 
ical point on transition mi = +1/2 <-» +3/2 at site la with 
n delays are labeled on the fit line 



point transition. 

The biexponential decay is attributed to the existence 
of a "frozen core" 013. A frozen core is formed by the 
Pr ion detuning nearby Y ions from the bulk Y, thereby 
inhibiting mutual Y spin flips within the frozen core. A 
strong frozen core is expected due to the Pr magnetic mo- 
ment ratio being an order of magnitude larger than that 
of Y. Although a frozen core has not been reported in 
Pr 3+ : YoSiO ^, it has been observed in directly analagous 
systems |l3j. On short time scales the rapid spin flips of 
the bulk Y ions dominate the dephasing. Since the bulk 
Y are well removed from the Pr ion they only weakly per- 
turb the transition frequency. On longer time scales the 
frozen core Y spin flips dominate the dephasing as they 
induce larger but less frequent frequency shifts. Given 
that the slow spectral diffusion has not plateaued for the 
longest T2 used, the correlation time of the frozen core is 
comparable or longer than 0.5s. 

Further investigation of the three pulse echos on the 
critical point transition (figure^J shows that for all values 
of n only a single, slow exponential decay was observed. 
Therefore dephasing of the critical point transition dom- 
inated, at all time scales by the frozen core Y spin flips. 
The fact that the bulk Y dephasing contribution is ab- 
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sent for the site la mj = +1/2 <-> +3/2 transition and 
present for the other two transitions is consistent with 
the dephasing at the critical point being due to second 
order Zeeman interactions. This suggests that the ap- 
plied field is accurate to the order of a Pr-Y interaction 
strength from the ideal critical point. Further optimisa- 
tion of the applied field is therefore unlikely to produce 
a significant gain in coherence time. 

Using the method described above a regime has 
reached where the correlation time of the dephasing in- 
teraction is extremely long compared to the duration of 
readily achievable RF pulses. This is exactly the regime 
required for effective implementation of quantum error 
correction schemes [TEl ITfil Il7j . The simplest implemen- 
tation shown to be sufficient for error correction in solid 



state qbitsll5( consists of a series of hard rephasing pulses 
(it) separated by the cycling time. In the present case the 
difference between a transition specific hard RF ir pulse 
(> 1/is) and the Tm would allow an effective error cor- 
rection scheme to be applied with a duty cycle < 0.1%. 

We have shown that a phase memory time of 82ms can 
be achieved for Pr 3+ ^SiO*, by using an external mag- 
netic field to minimise the transition sensitivity to mag- 
netic field fluctuations. Directions to further reduce the 
residual dephasing mechanism have been identified. We 
have shown that even in hosts containing nuclear spins it 
is possible to obtain spin based qbits that have coherence 
properties suitable for sophisticated quantum computing 
demonstrations . 
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